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c/3 I Abstract 

■ We investigate collective excitations of an atomic Bose-Einstein con- 
^ i densate in a recently discovered regime [D. O'Dell et al., Phys. Rev. 

I Lett. 84, 5687 (2000)] of a balance between electromagnetically in- 

■ duced 1/r attraction and a short-range interatomic repulsion. The cor- 
^ , responding frequencies for monopole, dipole, and quadrupole modes are 

I— 1 1 calculated numerically in the zero temperature case. 

^ : PACS numbers: 03.75.Fi, 34.20. Cf, 34.80.Qb. 

a^ ■ 
o : 

CN . Theoretical as well as experimental studies of dilute atomic Bose - Ein- 

^ I stein condensates (BECs) have made an impressive progress since 1995, when 

O I quantum degeneration in alkali metal wapours was reached at the first time 

^ ■ [H . These systems are interesting especially due to a presence of an interaction 

^ • between atoms weak enough to provide a development of quite precise theories 

^ , describing the properties of these non-ideal degenerate Bose-gases. The inter- 

^ I atomic potentials act on a short range, moreover, at ultralow temperatures 

C3 I they are effectively described by contact pseudo-potentials containing a single 

O ■ parameter, the s-wave scattering length a 0. Here we consider a > 0, i.e. the 

^ ■ short-range potential provides mutual repulsion of atoms. 

• ^ , Recently, a new kind of atomic BEC has been proposed in Ref.p|. Namely, 

^ i it is shown that the particular configurqation of intense off-resonant laser 

d I beams gives rise to an effective 1/r interatomic attraction between atoms lo- 
cated well within the laser wavelength (since inelastic collisions destroy very 
rapidly BECs of number densities more than 10^^ cm~^, a long- wavelength in- 
frared laser has to be considered). The directions and polarizations of different 
incident laser beams can be chosen in such a way that the static (near-zone) 
1/r^ component of dipole-dipole interaction energy is effectively "averaged 
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out", and only 1/r attraction remains. Finally, the two-particle interaction 
energy including both usual short-range interactions and such an " electromag- 
netically induced gravitation" reads as 

^i2(r) = -J^^i^) - (1) 

where M is the atomic mass and the long-range interaction constant u is 
expressed in SI units as B 

niqW 

The laser beam intensity and wavenumber are denoted by / and q, respec- 
tively; a is the atomic polarizability at the frequency q/c (practically, because 
the external radiation belongs to the IR range, a is equal to the static polaqriz- 
ability). A useful parameter = fi^ /{Mu) can be called the Bohr radius for 
the gravitation- like potential. Practically, the "gravitation" is weak enough, 
so always a^.^ a. 

In Ref. it is shown that a self-bound solution of the modified Gross - 
Pitaevskii equation taking into account the gravitation-like interatomic poten- 
tial exists. A cloud of atoms located in a radiation field of necessary geometry 
giving rise to 1/r attraction is stable without any confining external potential 
of magnetic or optical trap. The situation is similar to that of astrophysics 
where we have stars — objects those are stable due to a balance between 
gravitation and a Boltzmannian gas pressure. In our case the pressure arises 
from the two sources: interatomic repulsion due to the short-range potential 
and the "quantum pressure" corresponding to the kinetic energy term in the 
Schrodinger equation. 

Let us consider a zero-temperature case. Expressing the EEC order param- 
eter through condensate density n and phase (f) as usual [@], \l/ = nexp(i(/)), we 
rewrite the complex Gross - Pitaevskii equation as a set of two real equations: 

dn „ , ^ 

_ + V(nv) = 0, (2) 

where the local "gravitation potential" is 

U,{r) = -u f dh' 

J |r — r I 

The velocity field is introduced by the relation v = hVir/M. 



If ^ \Jaj,/a {N is the total number of atoms in the BEG) then the set of 
Eqs.(ll) admits an analytical steady state solution 0. Indeed, in such a limiting 
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case the kinetic energy term is small and can be omitted; the "gravitation" is 
balanced by the s-wave scattering. The density of such a self-bound degenerate 
bosonic cloud is equal to 



sin(7rr/i?o) 
4R| r 



eiRo 



0, 



(3) 



where 9 is the Heavyside function, _Ro = 'n'y/a^ is the radius of the cloud 
boundary (in Ref. there is an arythmetic error in the numerical coefficient 
in the expression for Rq). 

In the present paper, we calculate the eigenfrequencies of oscillations of such 
a self-bound cloud. We introduce small perturbations ni and vi of the density 
and velocity, respectively: n = Uq + ni, v = Vi. procedure of linearizing of 
Eqs.(|D is standard 0]. We choose the time dependence of perturbation in the 
form of exp{—iujt) . We calculate the eigenfrequenciues u within the Thomas - 
Fermi approximation, i.e., we neglect the kinetic energy term in the equation 
for the velocity. After elimination of the variable vi, we get finally 



u;V(r) + V 



smiri ~ / 

Ui r 
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(4) 



Eq.p) is written using dimensionless quantities 



TT 



U = U\ 



(5) 



It is convenient to introduce the "gravitation potential perturbation" variable 



f — r' I ' 



(6) 



so that rii = — V^^. Then Eq.(H) can be identically transformed to 



^2y2^ + V 



sm r 



0. 



(7) 



To solve this equation, firstly, we note that the orbital momentum / and 
its projection m to an arbitrary axis are conserved. So we can set 

where Yi^ is the spherical harmonic. The only remaining thing now is to set 
the boundary condition for Eq.(^. They can be derived from the tailoring of 
the inner and outer solutions at f = vr (i.e., r = Rq). First of all, a perturbation 
does not change the total number of atoms. So, for monopole modes, the 
radial derivative of the potential perturbation vanishes at f > vr. For modes 
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with / = 1,2,3, ... at the same distances ^/ decreases as r Hence, the 

boundary conditions are 

df 



— I- = 

\r=n ^ 



for the monopole mode and 



dr TT 



for the other modes {I = 1,2,3, ...). 

We solved Eq.(|^) numerically by expansion in series: 

^i(^) =J2'^ikjl{llkf/TT). 

k 

Here ji is the spherical Bessel function p. To satisfy the boundary conditions, 
one must choose the coefficient •jik equal to the k-th root of the equation 

jH7) = 0, (10) 

where 

/' = 1, / = 0, 

Z' = /-l, / = 1,2,3,.... 

For 1 = 1, the lowest frequency is exactly zero. Indeed, the lowest dipole mode 
corresponds to a translation of a self-bound BEC as a whole in a space free of 
any external trap potential. 

We calculated numerically the four lowest frequencies for / = 0, 1, and 2. 
The estimated numerical uncertainty is of order of one per cent. Thus, uj is 
equal to 

0.617, 1.27, 1.85, 2.39, ... for / = 0, 
0, 0.898, 1.51, 2.06, ... for / = 1, 
0.391, 1.10, 1.71, 2.28, ... for / = 2. 

Finally, we note that calculation of quantum depletion in a BEC self-bound 
by "electro magnetically induced gravitation" could be interesting. 
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